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Aircraft interference patterns in troposcatter signal records betray the plane's velocity 
across the link axis. The patterns also reveal the phase variations in the normal signal 
path. Similarly, the spectra of records free of aircraft patterns reveal the cross-axis velocity 
of winds aloft. The amplitude distributions of such records often come close to the Ray- 
leigh model but do not follow weather parameters closely. The distributions of ratios and 
products of correlated amplitudes also fit the Rayleigh model in records free of aircraft 
reflections. 



1. Introduction 

Interference between the normal signal and re- 
flections from aircraft flying in the common volume 
seen by transmitter and receiver leads to charac- 
teristic aircraft patterns in records of the signal 
level received over troposcatter links. These re- 
flections sometimes disturb TV reception, too. The 
interference patterns occur frequently in records 
made of continuous wave transmission at 10 kw and 
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near Boston to Schenectady, New York. 1 Usually 
the patterns seem a nuisance, but they may be worth 
a second look. They carry information on the air- 
craft^ velocity across the link axis. Moreover, 
they represent an extreme case of a single predomi- 
nant scatterer in the midst of scattering from many 
weaker sources, so they help to clarify features of the 
spectrum and the distribution of the signal ampli- 
tude. 

2. Aircraft Patterns 

To define our notation, we refer to a sideview of 
the link in figure 1. The aircraft at P, a distance 
T from the transmitter and R from the receiver, 
moves along a horizontal line p with velocity v and 
at a height h above the link axis. If p is projected 
down on to the horizontal plane through the link 
axis it will be found to cross it at a distance x from 
the axis midpoint, and at an angle \j/. The total 
signal path L by way of the aircraft is the sum T-\-R, 
and for its rate of change we find : 
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Figure 1. Notation for airplane patterns. 



under the assumption that h and p sin \p are much 
less than 1/2. The second time derivative of L is 
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Similar formulas are given by Pokorny in an article 
in Czechoslovakian not readily available [Pokorny, 
1961]. If the aircraft contributes a wave A and 
the normal path contributes a wave N to the receiv- 
ing antenna, as sketched in the lower part of figure 
1, the receiver output / represents the vector sum 
of N and A. We take N as the phase reference, 
and we use £ to indicate the phase difference between 
N and A. As £ changes, / traces the pattern char- 
acteristic of aircraft interference. If, for example, 
1-4 1 is less than \N\, as in figure 1, the pattern gives 
us the original wave amplitudes \N\ = 1 / 2 (I + J rI-) 
and \A\=){(I + —IJ), where I+, /_ are the extremes 
oil. 
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Figure 2. A typical airplane pattern and its analysis. 



In figure 2, one analysis of a sample 10 sec of 
interference has yielded the variations of \A\ and 
|iV| plotted beneath the record. In figure 3, the 
first graph includes these 10 sec, with A 2 plotted 
against time on a log-log scale; the second graph 
comes from a similar record. Time in each case is 
measured from the apparent center of the pattern. 
Both curves have knees — at about 14 and 6 sec 
respectively — and they fall off beyond the knees at 
slopes m of 6 and 10 respectively (ignoring the fine 
structure in the first pattern). Thus they can be 
approximated by the formula P^(l + [t/t ] 2 )~ m/2 . 

Now the Booker and Gordon [1950] formula for 
scattering from transparent inhomogeneities indi- 
cates that the power scattered will be proportional 




13 26 1.56 3.13 

SECONDS FROM CENTER 



6.25 



Figure 3. Power in the signal scattered by an airplane, No. 1 
and No. 2. 



to sin - ™ 0/2, where 9 is the angle at P in figure 1, 
between T extended and R, and m is about 4. For 
angles in the range of interest we have 



Psm~ m d/2~(l + [z/h] 2 )~ m/2 . 



(3) 



Here h is the height of the scatterer, and z=p sin \f/ 
is its horizontal distance from the link axis. If (3) 
is plotted on log-log scales for various values of m, 
we obtain figure 4, which consists of curves resem- 
bling those in figure 3. 
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Figure 4. A model of power scattered against off-axis position. 
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Reflections from shiny and opaque objects, how- 
ever, clearly will not depend on 6 in the simple manner 
of P in (3). The plane's velocity across the axis, v 
sin \f/, may be found from patterns like that in figure 
2, and then z k at the knee of the corresponding- power 
curve may be obtained. We find 4,700 and 3,500 ft 
off axis, respectively, for z k in the patterns analyzed 
in figure 3. These values of z k would equal the 
plane's height above the link axis, if (3) applied to 
scatter from aircraft; but the plane must, of course, 
be above the bottom of the common volume at 4,500 
ft (calculated from an earth with four-thirds the true 
radius). Again, these two values of z k fall very much 
short of indicating the half- width of the common 
volume: the receiver with an 18-ft dish and the 
transmitter with a 28-ft dish "see" regions 40,000 
and 26,000 ft in half-width respectively at midlink. 
Further data will be required to interpret correctly 
the shapes of the power curves in figure 3. 

We now obtain the velocity of the aircraft across 
the axis. Returning to figures 1 and 2, we note that 
the period T between adjacent minima of a pattern 
of interference can be fixed fairly accurately. If M 
is the length of the normal signal path, the difference 
(L — M) will change by a signal wavelength, X=0.328 
m, between adjacent minima: 



d(L-M) \ 
dt ~T 



(4) 



Except when the plane is directly above the link 
axis, or moving parallel to the axis, the chief term 
in (4) will be dL/dt rather than dM/dt, since the 
plane's velocity much exceeds that of the winds aloft. 
If we plot X/7 7 from a pattern such as appears in 
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randomly by the small effects of dM/dt. As shown 
in (2), the slope of this line represents the square of 
the cross-axis velocity of the plane. In figure 5, \/T 
from four records is smoothed out and plotted against 
the time measured from the pattern center. Typical 
aircraft velocities are obtained. 

Having identified the component dL/dt in plots of 
X/T, we turn next to dM/dt, and consider the phase 
stability of the normal path. It is difficult to moni- 
tor the phase of a tropospheric signal; neither refer- 
ence to very stable oscillators at both path terminals, 
nor the retransmission of the signal back to the 
transmitting site seemed feasible for us [Angell, Foot, 
Lucas, and Thompson, 1958]. We did, however, 
study the frequency dependence of the phase (or of 
the path length M) through data obtained by modu- 
lating the 915 Mc/s signal at 1 Mc/s. This gave 914 
and 916 Mc/s sidebands, hereafter labeled a and c, 
while b refers to the carrier frequency. The \/T 
analysis can be made of each component in records of 
aircraft interference; figure 6 gives an example. For 
the 14 sec analyzed from records of the simultaneous 
amplitudes of the a, b, and c signals (separated by 
filters), we find apparently unrelated variations of 
diVf(/, t)/c)t around a line of constant slope repre- 
senting d 2 L/dt 2 . 

Continuing the analysis of this record, figure 7 
gives in the dashed lines in the upper graph the 
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Figure 5. Analysis of four representative aircraft patterns. 



differences b(Mi—M b )/dt, i=a or c, between eacli side- 
band and the carrier b, for 7 of the 14 sec used in 
figure 6. Integrals of these curves from an arbitrary 
point in the —17th sec give the curves "upper" and 
"lower" in the lower graph of figure 7. In this 
record the frequency dependence of \d is small, i2 cm 
in 1^4 miles or two parts in 10 7 , for a change in 
wavelength of one part in 10 3 . The "sum" curves 
in figure 7 represent the second difference Q(J } t) and 
its integral: 
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The integral of Q may be recovered at the receiver 
without a "comparison" signal reflected from air- 
craft. The circuit used for this operation, designed 
by Dr. L. G. Abraham, Jr., mixed each sideband 
with the carrier to obtain two signals at 1 Mc/s, and 
it compared instants of zero-crossing of these 1 
Mc/s signals. The output was the second phase 
difference : 
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The phase </> a refers, of course, to the lower sideband 
written in the form/ a (£) cos (co a £+<£ G ), and similarly 
for ft and <j> c . A typical record of $>(f), about a 
minute long, made on August 12, 1959, appears as 
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Figure 6. Variations of X/T at three frequencies over the same fifteen second interval. 
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Figure 7. Phase fluctuation differences between sidebands and carrier. 



the top trace in figure 8. The other three traces in 
this figure give the amplitudes of signals a, b, and 
c. The arrows on these traces indicate the points of 
deep fade; the trace of <£> shows rapid changes at 
these instants. The discontinuities in the trace 



occur wherever <3> shifts over 180° either side of the 
(arbitrary) origin, and represent a periodicity in 
the instrument. 

The variation of <f> during this relatively calm 
period approached 720°. Changes of more than 
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Figure 8. Records of three envelopes, and second difference of phases. 



60=1 mm 



720° in a second occur in records of greater disturb- 
ance. Nevertheless, these records show that 
d 2 M/dt 2 changes gradually; sudden and wide fluctu- 
ations in $ like 180° in a millisecond would give 
spikes such as occur only when a signal vector 
passes near zero in a deep fade. 

3. Spectra 

t Much of the analysis and notation used to explain 
aircraft patterns helps also to describe the spectra 
of the amplitude records of signals not contaminated 
by aircraft reflections. The spectra we consider 
here were obtained by running several minutes of 
record, mixed with an audio signal, through a 
narrow filter (1 c/s wide, centered at 50 c/s) and 
integrating the filter output. Some records were 
processed for us at the University of Texas. These 
spectra appear to us to have resulted from scattering 
by a number of centers in fairly uniform motion 
across the path axis. The centers are presumably 
"blobs" of uneven dielectric constant in the atmos- 
phere, although their conformation must vary as 
much as the shapes of visible clouds, and may include 
rough or rippled layers as well as spheroids. The 
amplitude of radiation received from a blob will, 
of course, depend on its size and intensity, but also 
on the scattering angle, in accordance with (3). 
Thus for z/h larger than 1, the blob's contribution, 
beating against the average of a group of other 
wave vectors, will fall in amplitude as its zjh in- 
creases. The spectra as a whole tend to take the 
form of figure 4. Figure 9 shows examples, for a 
succession of 7 min recorded May 8, 1958, and 
reproduced here to show the consistent appearance 
of a knee, and a drop-off with exponent about 
m=4. The 8th min, contaminated by aircraft, 
has a different and typical shape, a plateau or 
bump on the drop-off. 

The "blobs" drift with the air movements, and 
scatter wavelets that differ slightly in frequency from 



the transmitted frequency f b . The Doppler shift 
depends, of course, on the rate of change of the 
path-length M via the blob; 

A/=/ _ /6 4^=i^M- 2 , (8) 
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assuming a horizontal drift and (2a://) 2 <^<^l. If the 
blob happens to move so M is constant, that is, on 
the surface of an ellipse of revolution with foci at the 
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Figure 9. Power density spectra from May 8 t 1958. 



080071— G3- 



409 



two antenna sites at the link ends, then A/=0, but 
when the blob moves perpendicular to such a surface, 
Af is a maximum for a given speed [Laaspere, 1958]. 
Usually, and particularly in winter, the winds carry 
scatterers horizontally, and we may solve for their 
cross-axis velocity, using for A/ the frequency at the 
knee of the spectrum, our link dimensions for IX, 
and for h=z the estimate of 4,000 m: 



dz/dt=v sin \l/=4A2 Af m/sec. 



(9) 



Figure 10 shows four spectra from records taken in 
December 1957, in extremes of weather. From 
Af at the knee of each curve we obtain cross-axis 
velocities of the scatterers. In figure 11 we compare 
these with actual velocities of winds aloft. The 
latter come from soundings [U.S. Department of 
Commerce, Weather Bureau, Daily Series, Synoptic 
Weather Maps, Part II, 1957] at the Albany Weather 
Station, obtained at 7 a.m. or 7 p.m. Although 
Albany is on the link axis, it lies 60 miles from the 
axis midpoint; moreover, our records were made 
between 8:30 a.m. and 4:30 p.m. so the Albany 
weather data are not ideal. Fortunately, winds 
aloft are fairly persistent phenomena, characteristic 
of large areas. In figure 11, points No. 1 through 
No. 4 correspond to curves so numbered in figure 10. 
The extreme point No. 1 falls nearly on the diagonal 
line of complete agreement between spectral and 
weather velocities. Points No. 7 and No. 10 also 
fall on it, while No. 14 is the worst case. When No. 
14 was recorded, March 26, 1958, the winds aloft, 
although between 30 and 40 m/sec, were almost 
parallel to the link axis. A change of 10° in wind 
direction would put No. 14 on the diagonal. 
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Figure 10. Spectra from extremes of weather. 
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Figure 11. Corr elation of cross winds and spectrum. 

Although for the aircraft patterns we found m 
from 6 to 10, most spectra drop off at a slope from 
4 to 6. The slope of the spectra must depend on the 
distribution of scatterers in height and on the 
antenna apertures if these are narrow, according to 
suggested modifications of the Booker-Gordon for- 
mula [Villars and Weisskopf, 1955]. 

4. Distributions 

The addition of a large number of wave vectors, 
none of them predominant, arriving with random 
phases from an assortment of scattering centers, 
yields a resultant wave vector whose amplitude tends, 
as a consequence of the central limit theorem [Lawson 
and Uhlenbeck, 1950], to have the Rayleigh dis- 
tribution 



P(r) dr=(r/a 2 ) exp (—r 2 /2a 2 ) dr. 



(10) 



An intensive analysis of 2 l / 2 hr of our records by R. G. 
Finney [1958] showed sections, particularly in records 
from calm days, with long periods of stationary rms 
levels, and distributions satisfying ax 2 test for a 
Rayleigh population, at the 5 percent level of proba- 
bility of rejecting the hypothesis, when the distribu- 
tion was, in fact, Rayleigh. 

We also collected samples of the cumulative distri- 
bution, obtained by reading the recorded signal 
envelope into a set of clocks controlled by relays 
triggered at preset levels (a Gates multilevel recorder). 
The results for seven consecutive minutes of clear 
record and 4 min (dashed curves) contaminated by 
aircraft are given in figure 12. There we plot the 
cumulative distribution, the fraction P of the record 
that lies below the amplitude r (abscissa), against X= 
r/r m (ordinate), where r m is the median signal level. 
Curve No. 1 coincides with the Rayleigh model, 
and the succeeding minutes yield curves not far 
from it, until the aircraft interference begins. 
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If we have interference between two equal scat- 
terers, separated a distance d, and at a height h, 
drifting at a velocity v across the link axis, the signal 
reception may take the following form : 

I(t)~(f+g 2 -2jg cos uf) 1 *; (11) 

where 

i//=i + (^/A) 2 , 1/0=1 + (vt-dy.hr\ 

and 

u=2tt/\. (d(M f -M g )/df) . 

The distribution of 1(f) for d/h=6, accumulated in 
the interval 0<C^<Cd, is plotted as the dotted curve 
in figure 12, and may help to explain the curve for 
the eleventh minute. 

Figure 13 shows distributions for the records from 
December 1957, from which the spectra in figure 10 
were taken. Clearly the distributions give no indi- 
cation of the weather: the 20th and 26th had strong 
winds aloft and high frequency knees in their spectra, 
while the 13th and 18th are calm days, yet the 18th 
and 20th have similar distributions and the 13th and 
26th are also paired. 

Departures from the Rayleigh distribution can 
sometimes be represented by a single scatterer pre- 
dominant among small random scatterers. The lines 
in figure 14, from a numerical intergration by Norton, 
Vogler, Mansfield, and Short [1955], show distribu- 
tions for various values of k, the ratio of the rms sum 
of small scattering vector amplitudes to the pre- 
dominant vector's amplitude. (The distribution is 
Gaussian on the wast- vector plane, but the center 
is displaced from the origin by the predominant 
vector.) Data from a record without obvious air- 
craft patterns, taken April 1958, appear as dashed 
lines on figure 14; the original data were close to 
Eayleigh above median amplitude, but below Ray- 
leigh for lower amplitudes. Adding 4.3 /iv (13% 
of the median 33 /xv) to each point moved the data 
to overlap the k=l curve over its whole length. 
This type of change enabled us to fit many curves 
to the Norton model. However, in each case the 
voltage increment is different, and does not, appar- 
ently, correct an error in data. We merely note 
that some distributions depart from modified Rayleigh 
as if by displacement of the zero level. 

In studies of the bandwidth of the tropospheric 
path, a parameter of interest is the correlation p 
between received amplitudes of signals separated 
slightly in frequency. This correlation is subject to 
statistical fluctuations, and depends also on weather 
conditions in the common volume. One can process 
the amplitude records to give the average values, 
<r!r 2 > or ^rjr^^, and with these determine p 
from formulas such as we give in an appendix ; alter- 
natively one can find the cumulative distributions of 
?\r 2 and ri/r 2 , and determine the correlation by com- 
parison with curves plotted in figure 15. An example 
of the second method, which provides a consistency 
check in the shape of the curve, appears in figures 16 
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Figure 12. Amplitude distributions for eleven minutes, May 
8, 1958. 
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Figure 13. Amplitude distributions in extremes of weather. 
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Figure 14. Distribution of amplitude of sum of a constant 
vector and random vectors. 



and 17. The data analyzed in this way from records 
made August 5, 1959, are typical of a dozen so proc- 
essed. They represent here about 4 min of normal 
record, and 1 min of aircraft interference. 

The normal 2 distribution at all three frequencies, 
a, b, and c, is close to Rayleigh in figure 16, and the 

2 "Normal" is used here as "commonly occurring," rather than normal in the 
statistical sense. 
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Figure 15. Distributions of functions of two correlated Rayleigh amplitudes. 



ratios are distributed along curves parallel to the 
models in figure 15. 3 In figure 17, the products from 
the normal records also follow the model curves, but 
the products involving aircraft patterns, as might 
be expected, have more and higher peaks. However, 
again in figure 17, the ratios from aircraft patterns 
are less disturbed than the products. Finally we 
compare the correlation parameter R, determined 
from <Crir2>, with p 2 , determined from the cumula- 
tive distribution of rjr 2) for each pairing of the three 
normal records : 



Pair 


B 


P 2 


(a,b) 


0.32 


0.30 


(b,c) 


0.49 


0.50 


(c, a) 


0.25 


0.20 



(The definition of JS, and its relation to p are given in 
the appendix.) 



3 The ratio distributions are "double": points for X less than 1 are plotted as 
1/X, and the sign of is changed, so the distributions of T\\ri and r 2 /n are not dis- 
tinguished. In general, the data are symmetric about the median, so that points 
XSl and mirror points X<1 fall on the same curve, and are averaged in plotting 
figures 16 and 17. 



5. Summary 

Airplane interference patterns in records of signal 
amplitude betray the plane's cross-axis velocity, and 
the phase variations of the normal path. The inter- 
ference is an exaggerated form of scattering from a 
predominant scatterer in the presence of random 
signals from many small ones. It markedly affects 
both the spectrum and distribution of the envelope 
record. The spectrum in the absence of aircraft can 
be related to a model of scatterers drifting with the 
winds aloft. The correlation of signal amplitudes in 
nearby frequencies fits a Rayleigh model for the joint 
distribution, as represented by derived distributions 
for the product and ratio of signal amplitudes. 



The receiving equipment and the phase comparator 
were designed and assembled under the supervision 
of Dr. L. G. Abraham, Jr., who also participated in 
the studies of correlation and effective bandwidth 
and other parts of the project. J. S. Brookman and 
T. E. Kotary processed the data. 



412 



4 — 



LU 
Q 

Z> 

< 



o 



i — r 




RATIOS 
NORMAL RECORD 



I 



i r 



-RAYLEIGH NORMAL RECORD 



V\ ^" A 




0.3 



0.95 0.9 0.8 0.7 0.95 0.9 0.7 

P= FRACTION OF TIME ORDINATE NOT EXCEEDED 
Figure 16. Distributions of amplitudes and amplitude ratios at three frequencies, taken August 5, 1959. 
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Figure 17. Ratio and product distributions from August 5, 1959. 



7. Appendix 

If two distributions derive from Gaussian popu- 
lations on the complex wave vector plane, the cor- 
responding joint distribution can be written with 
correlation parameters p and ^ [Wheelon, 1959], 

^ 4 rfl2=exp —-4 rl+r 2 2 — 2r 1 r 2 p cos (0i— 2 +tfO V 

1 /y> ryt 

i 2 Xti« —-^i • dr x dr 2(10462. ( A. 1 ) 

1 — p _J 47T 1 — p 

Here the complex wave vectors, x i -\-jy i =r i exp jd t (i= 
1 or 2), are normalized so that <£?> = <J/?> = 1, 
and p cos t=<Cx l x 2 >=<yiy2>,f sin ^=<ai2/ 2 >= 
— < C.x 2 yC> 1 the brackets indicating statistical aver- 
ages. Integration over 6 Y and over 2 leaves the 
joint distribution of two Rayleigh-distributed ampli- 
tudes: 



p%{r u r 2 ; pyir^ 

-^[-K©)Kt^i)^ «* 

Now pAtr u r 2 ,d 1 ,d 2 ; p, ^) in (A.l) is periodic in (0 A — 2 ). 
But in our case the difference in path lengths, 
represented by 0i — 2 for the two signals, will be 
distributed to peak at some value near zero, rather 
than to show periodicity. (A Gaussian distribution 
for the phase difference might be a better representa- 
tion than £> 4 ; an attempt, not wholly successful, to 
derive such a representation appears on page 512 
of Hirai, Fukushima, and Kurihara [I960].) This 
objection, however, does not apply to the derived 
expression, p 2 , which has no phase dependence; 
moreover, it reduces to the same Rayleigh form on 
integration over r x or r 2 . It appears to provide an 
adequate model for much of our data on correlated 
Rayleigh amplitudes. 
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From (A. 2), after integration over the product 
/Vj in the interval (0, <») and over the ratio r x /r 2 
in the interval (1, X), we obtain: 



0(p)= 



sinh a 



[cosh 2 a— p 2 ] 



211/2 



> with 1 + /3 

=2P,expa=X (A.3) 



From (A.2), after integration over ri/r 2 in (0,co) and 
over Tir 2 in (0, ?/), we obtain: 

P(y) = j\ 1 ~P 2 )K (?/) 7o(p2/) ||^, 

= l-y[I (py)K l (y)+pK (y)I 1 (py)}; 



nr 2 

"1-/ 

(A.4) 



In figure 15 we plot the expressions P(X) from (A.3) 
and (A.4) for four values of the correlation p. In 
plotting (A.4), the ordinate X has been normalized 
(dividing y by its median v alue) , so P ( 1 ) = % . 

From (A.2) we can also determine the expected 
values: <nr 2 >=2E( P )-(l- p 2 )K(p), and <r } /r 2 > 
=E(p). Here E and K are complete elliptic in- 
tegrals of the first and second kinds (while 1 , I\, 



K and K\ in (A.4) are Bessel functions). The ex- 
pressions in question are conveniently modified to 
emphasize their dependence on p. We are led to 
define correlation parameters R and S thus: 



R- 



<nr 2 >-«r» 2 
- <r 2 > _ {<r>) 2 



=0.91()p 2 (l + p 2 /U)+.. .), 



and 



S-- 



<r 1 /r 2 >-<r><l/r> 
!-<>><l/r> 



-0.688p 2 (l+3p 2 /16+ . . .)• (A.5) 



One expects p to be a monotonic decreasing func- 
tion of the frequency difference between signals, 
so p(a, c) should be less than p(a, b) or p(b, c). How- 
ever, there is apparently no necessary relation be- 
tween these three correlations. Moreover, in a 
search for closed expressions for the integrals in- 
volved in <Cr a r b r c ^>, and in the cumulative dis- 
tribution of the triple product, P(r a r b r c ), we have had 
little success. 

(Paper 67D4-276 
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